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Production of high-energy /i neutrinos from young neutron stars 
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Young, rapidly rotating neutron stars could accelerate protons to energies of ~ 1 PeV close to the stellar 
surface, which scatter with x-rays from the stellar surface through the A resonance and produce pions. The pions 
subsequently decay to produce muon neutrinos. We find that the energy spectrum of muon neutrinos consists of 
a sharp rise at ~ 50 TeV, corresponding to the onset of the resonance, above which the flux drops as up to 
an upper-energy cut-off that is determined by either kinematics or by the maximum energy to which protons are 



accelerated. We predict event rates as high as 10-50 km 
fluxes would be detectable by IceCube. 



1. Introduction 

A new window in high-energy astronomy is 
opening as existing neutrino detectors are im- 
proved and new ones are developed. Since neutri- 
nos produced in astrophysical systems are unim- 
peded by interstellar matter on their way to 
Earth, detections will provide a new way to study 
the highest energy phenomena in the Universe. 
Astrophysical neutrinos are expected to arise in 
many environments in which neutrinos are pro- 
duced by the decay of pions created through 
hadronic interactions (pp) or photomeson produc- 
tion (jrf). Neutrinos may be produced by cosmic 
accelerators, like those in supernova remnants 
active galactic nuclei micro-quasars |Sj and 
gamma-ray bursts |4l5j . To detect these neu- 
trinos, several projects are underway to develop 
large-scale neutrino detectors under water or ice. 
AMANDA- II 0, in the South Pole, and Baikal 
[7| are the two neutrino telescopes currently run- 
ning, whereas ANTARES [HI and NESTOR ^ 
are under construction in the Mediterranean Sea. 
Those telescopes belong to a first generation with 
instrumented volume smaller than 0.02 km^. The 
IceCube detector ^UIj with a volume of about 
1 km'^, is under construction on the same site 



from relatively young, close neutron stars. Such 
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of AMANDA. The NEMO project |TT| is in its 
starting phase, and will be a cubic kilometer size 
detector located at Capo Passero, Southern Italy. 

As neutrino astronomy comes of age, it is im- 
portant to get some idea of what the sources 
might look like to aid in their detection. Recently, 
we proposed |12I13) that young {tage ^ lO^yr) 
and rapidly-rotating neutron stars could be in- 
tense neutrino sources. Here we summarize this 
work. 

2. The Model 

Neutron stars have enormous magnetic fields 
(> 10^^ G, typically) and high rotation rates 
(tens of Hertz), causing them to act as very 
powerful unipolar generators. Charges stripped 
off the highly-conductive surface are accelerated 
somewhere above the stellar surface. As the 
charges are accelerated along magnetic field lines, 
they emit curvature radiation, which scatters 
with the magnetic field and produces an e^e~ 
cascade. The cascade produces a radio beam 
along the magnetic field lines that open to in- 
finity. If the beam angle and direction to Earth 
are favorable, the star can be detected as a radio 
pulsar, with one radio pulse per rotation period 
(in some cases, two pulses, when emission from 
both magnetic poles is seen). 
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If the stellar magnetic moment has a compo- 
nent anti-parallel to the spin axis (half of neutron 
stars) , ions will be accelerated off the surface (see 
Fig.l). If energies of ~ 1 PeV per proton are 
attained, pions will be produced through pho- 
tomeson production as the protons scatter with 
surface x-rays, producing a beam of /i neutrinos 
with energies above 50 TeV. Detection of such 
neutrinos would be a fascinating discovery in its 
own right, and would provide an invaluable probe 
of the physical conditions that prevail in the mag- 
netosphere of a neutron star. 

For photomeson production to occur, ions must 
be accelerated to very high energies close to the 
stellar surface. How large might the accelerat- 
ing potential be? Goldreich and Julian developed 
the first model of a quasi-static magnetosphere 
|14| . By assuming a dipolar configuration with 
magnetic axis parallel to the rotation axis, they 
showed that the potential drop across the field 
lines of a pulsar with angular velocity f2 = 2tt/p 
(where p is the period) from the magnetic pole 
to the last field line that opens to infinity is of 
magnitude 

^ ^ 7 X 10i«i?i2i?6V„^, Volts. (1) 

Here B = 10^^i?i2G is the strength of the dipole 
component of the field at the magnetic poles, 
R = 10^i?6 cm is the stellar radius and p„^s 
is the spin period in milliseconds. In equilib- 
rium (not realized in a pulsar), a co-rotating 
magnetosphere would exist in the regions above 
the star in which magnetic field lines close; the 
charge density would be pq ~ eZno ~ B/pc 
(cgs units), where hq is the equilibrium number 
density of ions that would short out the compo- 
nent of the electric field along the magnetic field. 
Deviation from corotation will lead to charge- 
depleted gaps somewhere above the stellar sur- 
face, through which charges will be accelerated to 
relativistic energies |15llti| . Suppose that charge 
depletion occurs near the stellar surface, with a 
characteristic density fdno, where < 1 is an 
unknown depletion factor. If the neutron star is 
young, its surface will emit in soft x-rays, and the 
protons in accelerated nuclei will scatter with this 
radiation field. If the protons are sufficiently ener- 



getic, they will exceed the threshold for photome- 
son production through the resonance. The 
A+ quickly decays to a 7r+, and muon neutrinos 
are produced through the following channel: 

P7 A+ n7r+ — > nv^e^VeV^. (2) 

The proton energy threshold for A+ produc- 
tion is given by 

tpe^ > 0.3 GeVVg, /g = (1 - cose)-\ (3) 

where e-y is the photon energy and 6 is the in- 
cidence angle between the proton and the pho- 
ton in the lab frame. Young neutron stars typ- 
ically have temperatures of Too — 0.1 keV, and 
photon energies — 2.8kToo{l + Zg) ~ 0.4 keV, 
where Zg ~ 0.4 is the gravitational red-shift and 
Too is the surface temperature measured at in- 
finity. The proton threshold energy for the A 
resonance is then ep^th — ^Vikcv/ff PeV, where 
To.ikcv = (A:Too/0.1 keV). Therefore, if the po- 
tential along field lines is only 1% of the full 
potential A<i> across field lines in the equilibrium 
magnetosphere, protons can reach the A reso- 
nance. 

Before continuing, we mention that our as- 
sumption that PeV protons can be produced in 
the stellar magnetosphere is a strong one; it is 
generally thought that pair production near the 
stellar surface will quench the field and limit the 
potential drop along the field to 1 TeV jl7ll8| . 
These calculations are subject to numerous uncer- 
tainties, and apply only in equilibrium, not the 
case in a neutron star magnetosphere. Though 
potential drops as high as ~ 1 PeV might not be 
attainable, we consider the possibility sufficiently 
interesting to explore further and to possibly con- 
firm or refute with future neutrino observations. 
Even null results would provide a valuable probe 
of the poorly-understood neutron star magneto- 
sphere, confirming a limit that so far has only a 
theoretical basis. 

Assuming that a potential of order is avail- 
able for acceleration along field lines, a necessary 
condition for the A resonance to be reached is 

BuP^nlTo.ikcV > 3 X 10-4. (4) 

Taking To.ikcV — 1; typical of pulsars younger 
than 10^ yr, there are 10 known pulsars within 
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a distance of 8 kpc that satisfy this condition 
|19|. about half of which should have positively- 
charged magnetic poles; these are potentially de- 
tectable sources of fi neutrinos. The best can- 
didates are young neutron stars, which are usu- 
ally rapidly spinning and hot. In the photomeson 
production process of eq.(2), the muon neutrinos 
receive 5% of the energy of the proton. Typi- 
cal proton energies required to reach resonance 
are ~ 1 PeV, so the expected /i neutrino energies 
will be ~ 50 TeV. The conversion probability for 
a proton is small (< 1%); only a small number 
of protons are excited to the A resonance and 
the pulsar beam is essentially unaffected. More- 
over, since the accelerated protons are far more 
energetic than the radiation field with which they 
interact, any pions produced through the A res- 
onance, and hence, any muon neutrinos, will be 
moving in nearly the same direction as the proton 
was when it was converted. The radio and neu- 
trino beams will be approximately coincident, so 
that some radio pulsars might also be detected as 
neutrino sources. We see the radio beam for only 
a fraction /b of the pulse period, the duty cycle. 
Typically, /h ~ 0.1 — 0.3 for younger pulsars. We 
take the duty cycle of the neutrino beam to be 
fb (but see below). In Ref. [T3|, we estimated 
the phase-averaged neutrino flux at Earth result- 
ing from the acceleration of positive ions, at a 
distance d from the source, to be 

(I), - 2chfd{l - fd)no (^^^ Pc, (5) 

where fdXio is the space charge density in the de- 
pleted region, where proton acceleration is taking 
place, and Pc is the conversion probability. For 
the purposes of calculating the spectrum, we use 
the following differential neutrino flux: 

We are interested in obtaining upper limits on 
the flux and so take fd — 1/2 and Z = 1. The 
detailed calculation of the spectrum dPc/de^, has 
been reported in and we summarize these 
results next. 




Figure 1. For a neutron star with a compo- 
nent of its magnetic moment antiparallel to its 
spin vector (these two vectors are shown here as 
completely antiparallel for simplicity), the sur- 
face of the star will be charged as shown. In 
the magnetic polar caps (one shown here), pos- 
itive charges will be stripped off the surface and 
accelerated through a region of some dimension 
L (shaded region). We postulate that the charge 
depletion in this region is sufhcient that accel- 
erating voltages of ~ 1 PcV per proton exist, at 
least sporadically. The protons scatter with x-ray 
photons from the stellar surface, producing pions 
through the A resonance. 



3. Neutrino Spectrum 

In calculating the neutrino spectrum, we have 
accounted for the fact that the proton accelera- 
tion takes place over a finite distance above the 
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Table 1 

Estimated upper limits on the /i fluxes at Earth. Numbers followed by question marks indicate 
guesses. Radio pulsar spin parameters were taken from the catalogue of the Parkes Radiopulsar Sur- 
vey (www.atnf.csiro.au/research/pulsar/psrcat/). Temperatures and limits on temperatures were taken 
from the references indicated. The temperature upper limits on the Crab and J0205-I-64 were used. The 
integrated conversion probability Pc is reported for the case of linear acceleration, and assuming L = 0.1. 



Source 


^kpc 


age 
yr 




Bl2 


To.lkoV 


fb 


Pc 








dN/dAdt 

km^2 yr^ 


Crab 


2 


10^ 


33 


3.8 


< 1.7 [201 


0.14 


1.6 


X 


10 


-3 


45 


Vela 


0.29 


104-2 


89 


3.4 


0.6 EI] 


0.04 


7.2 


X 


10 


-5 


25 


J0205+64 


3.2 


102.9 


65 


3.8 


< 0.9 [22 


0.05 


2.4 


X 


10- 


-4 


1 


B1509-58 


4.4 


103.2 


151 


15 


1? 


0.26 


3.4 


X 


10- 


-4 


5 


B 1706-44 


1.8 


104.3 


102 


3.1 


1? 


0.13 


3.4 


X 


10- 


-4 


5 


B1823-13 


4.1 


104.3 


101 


2.8 


1? 


0.34 


3.4 


X 


10- 


-4 


2 


Cass A 


3.5 


300 


10? 


1? 


4 [22 


0.1? 


2.1 


X 


10 


-2 


50 


SN 1987a 


50 


17 


1? 


1? 


4? 


0.1? 


2.1 


X 


10 


-2 
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surface. We have considered an acceleration law 
of the form 

"-'fM <') 

where eg = 0.3 GeV^, z is the height above the 
stellar surface and L is the length scale over which 
the proton is accelerated before reaching sufficient 
energy to undergo the A resonance. We consider 
linear acceleration (7 — 1), corresponding to a 
constant electric field in the charge-depleted zone, 
and quadratic acceleration (7 — 2), correspond- 
ing to an accelerating field that grows linearly 
with height. We have also accounted for the fi- 
nite width of the A resonance in the cross section. 

The spin parameters of 8 potential neutrino 
sources are given in Table 1. In Fig. [21 we show 
the neutrino energy flux for the Crab pulsar, as- 
suming its surface temperature is equal to the 
observational upper limit. The Crab pulsar is lo- 
cated in the northern hemisphere and so is mon- 
itored by AMANDA-II. The spectrum calculated 
in our model is represented by the solid line, and 
it has been calculated by assuming linear accel- 



eration. For comparison, we also show two inde- 
pendent calculations obtained by Bednarek and 
Protheroe [23| (dotted line), and Guetta and Am- 
ato [21] (dashed line) in the case that neutrinos 
are emitted from the nebula. Those models pre- 
dict a low muon event rate in comparison to our 
model, as shown in Ref. |26| . In our model, the 
spectrum begins sharply at 

~ 7QT^^\ GeV, (8) 

corresponding to the onset of the resonance. 
At higher energies, the spectrum drops approx- 
imately as e~2, as the phase space for conversion 
becomes restricted; higher energy neutrinos are 
produced by protons that have been accelerated 
to greater heights, where the photon density is 
lower and solid angle subtended by the star (as 
seen by the proton) is smaller. At some maximum 
energy, the spectrum is suddenly truncated by ei- 
ther kinematics (solid curve) or the termination 
of the proton acceleration as limited by the mag- 
nitude of the acceleration gap (not shown, since 
this cut-off cannot be predicted). 
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Figure 2. The predicted neutrino energy flux for the Crab pulsar for hnear acceleration of the protons 
(solid line). Calculations from other models are also shown for comparison. 



4. Estimated Count Rates and Discussion 

Large-area neutrino detectors use the Earth as 
a medium for conversion of a muon neutrino to a 
muon, which then produces Cerenkov light in the 
detector. The conversion probability in the Earth 
is P^^-.;, ~ 1.3 X 10-S(e^/1 TeV), where is the 
energy of the incident muon neutrino "27". The 
muon event rate is 

^-fde^-^P (9) 

Estimated conversion probabilities for the Crab, 
Vela and 6 other pulsars are given in Table 1 (col- 
umn 8) for a characteristic acceleration length 
L = 0.1. The final column gives the estimated 
event rates. Complete consideration of the kine- 
matics gives event rates are a factor of 10 — 30 
lower than estimated in our previous paper |12j. 

Neutrinos are produced at relatively high rates 
only if the protons are accelerated through the 
resonance close to the star. In this case, we ob- 
tain integrated count rates of several to ^ 50 
km^^ yr^^. Such count rates should be easily 



detected by IceCube, and possibly by AMANDA- 
II or ANTARES with integration times of about 
a decade (IceCube is planned to have replaced 
AMANDA-II by then) for depletion factors of 
fd — 1/2. While the characteristics of the spec- 
trum presented here are robust, we caution that 
the events rates we obtain are very rough upper 
limits, subject to many uncertainties. For ex- 
ample, we have assumed that the neutrinos are 
beamed into the same solid angle as the radio 
beam, which might not be a correct assumption. 
The radio beam is thought to be produced at 
about lOi? 12^. In our model, the pions are pro- 
duced much closer to the star. They then prop- 
agate to ~ lOOOi? before decaying to neutrinos. 
At this distance from the star, the field is unlikely 
to be dipolar, and it is difficult to say anything 
definite about the distribution of pion velocities 
in this region. If the neutrinos form a beam, the 
beam may be more or less collimated than the 
radio beam. If it is more collimated, the neutrino 
event rates could be higher than estimated here. 
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Results of 807 d of data from AMANDA-II are 
now available UHl- AMANDA-II has detected 10 
events (over a background of 5.4) from the di- 
rection of the Crab pulsar, with energies higher 
than 10 GeV. This result, though intriguing, is 
not statistically significant; IceCube will be able 
to confirm or refute this result. While it would 
be more exciting to see neutrinos from pulsars, 
the accumulation of null results over the next 
decade would be interesting as well; it would 
probably mean that photomeson production is in- 
effective or non-existent in the neutron star mag- 
netosphere, thus providing a bound on the accel- 
erating potential that exists near the neutron star 
surface. 
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